How long will a bubble be? 
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A soap bubble is a metastable object that eventually breaks. Indeed, the soapy water film thins 
until rupture, due to drainage and evaporation. In our experimental investigations, floating bubbles 
at the surface of a liquid bath have been considered. Their lifetime has been measured and reported 
with respect to their radius. Large bubbles last longer than small ones. Moreover, small bubbles have 
more predictable lifetimes than large ones. We propose a general equation for that lifetime, based 
on the lubrication theory. The evaporation is shown to be an essential process which determines 
the bubble lifetime. 

PACS numbers: 47.15.gm, 47.55.D-, 47.55. dd, 47.57.Bc 



A child succeeds in blowing a soapy bubble, spherical, 
light, fragile. The bubble flies through the air, avoids 
some mortal hydrophobic obstacles and eventually dies 
after a frontal collision with a wild red poppy. How long 
would it have lived ? A second question should also be 
addressed: is the bubble lifetime related to its size when 
prevented from any accident ? Some indications can be 
found in the literature. The lifetime of large bubbles 
made of PDMS oil with a large viscosity (lOOOPa-s) have 
been studied in [l[ . This liquid had been chosen to avoid 
dust contamination and evaporation. The bubble life- 
time has been shown to be related to the thinning of the 
film at the top of the bubble. The thickness is found to 
decrease according to an exponential law until a critical 
value for which the film breaks. This allows to define a 
lifetime r for the bubble. In the case of PDMS bubbles, 
the lifetime has been found to be inversely proportional 
to the bubble radius. Although the thinning behavior 
has also been observed in surfactant-water systems Q, 
flows in soapy films are very different from the ones in 
the PDMS system. Firstly, surfactant molecules rigidify 
the interface: The flow cannot be described as a plug 
flow. Secondly, the evaporation cannot be neglected. In- 
deed, let us remind that when a soap film is prevented 
from any evaporation and nucleation by dusts, the film 
is metastable and the bubble may last as long as the 
experimentalist's patience. On the other hand, physical 
processes as drainage, nucleation of holes by dusts and 
evaporation reduce the lifetime of the bubble. Are those 
processes influenced by the size of the bubble? Some 
experimental works can be found about the lifetime of 
n-butanol and n-nonanol bubbles [3)I3]- In these papers, 
the motion of the bubbles through the liquid is shown to 
modify their lifetime because the surfactant molecules lo- 
cated at the interfaces are redistributed during the travel. 

In this letter, the lifetime of a single bubble of soapy 
water is investigated. The original part of this work is 
that evaporation plays here a determinant role in the 
bubble lifetime. In order to stabilize the bubble, surfac- 
tant Triton X-100 is added to water with a concentration 
of 2mMol/L. This corresponds to approximately 10 times 




FIG. 1: Pictures of bubbles at liquid-air interfaces. The up- 
per part of the bubble (above the bath interface) can be fit- 
ted by a sphere of radius R (Here approximately drawn) . (a) 
Bubble with R = 1.67mm (0.005mL in volume); the air bub- 
ble remains below the liquid-air interface, (b) Bubble with 
R = 31.6mm (80mL in volume); the bubble shape is an 
emerged hemisphere. 



the critical micellar concentration. Indeed, the lifetime of 
the bubble increases with the surfactant concentration, 
and even saturates for large concentrations [3, 0] ■ There- 
fore, using a liquid with a concentration of surfactant 
10 times above the CMC prevents from any effects due 
to the depletion of surfactant molecules. Moreover, inter- 
faces can be considered as nearly rigid ones and the plug- 
flow component is negligible. That surfactant molecule 
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is non-ionic and is used as a standard in many physico- 
chemical experiments. 

Contrary to the previous works P, H, 01 j bubbles are 
created at the surface of a water bath by using a syringe 
filled with air. In so doing, the distribution of the surfac- 
tant molecules along the surface of the bubble is uniform 
and relatively homogenous. The volume of the bubble 
is tuned by using different syringes, from 2fiL to 40mL. 
The relative error made on the volume is estimated to be 
less than 5%, mainly due to the graduation reading. The 
atmospheric conditions are constant during the whole ex- 
periment. In particular, the temperature is 22°C and the 
relative humidity ( = 0.54. In order to avoid dust con- 
tamination and hazardous variations of the evaporation 
rate, the system is set up in a constant laminar flow of 
purified air with a velocity about O.Olm/s (Captair Flow 
- Erlab). Moreover, that ensures the atmospheric per- 
turbation (wind) to be constant and the same for all our 
experiments. A digital camera records the bubble from 
the side (Fig[T]). The relative error on the radius mea- 
surement is less than 2% (due to the resolution of the 
camera). Moreover, an average of the geometrical prop- 
erties is made over three different bubbles with the same 
volume. 




FIG. 2: Shape of a bubble near a liquid-air interface. Radii of 
curvature at the bottom Rb and at the top R are related to 
the bubble depth zb according to Eq.(IT}. The dotted curve 
corresponds to the contour along which pressure variations 
are integrated. 



The shape of a bubble (or a droplet) near a liquid-air 
interface is the result of a competition between reduced 
gravity (gravity minus buoyancy effects) and surface ten- 
sion effects 0. When gravity is negligible compared to 
surface forces, the bubble remains spherical and fully sub- 
merged below the perfectly planar surface of the bath 



(Fig 1(a)). This configuration minimizes the surface en- 
ergy of the system. When gravity is much more impor- 
tant than interfacial tension, the bubble is nearly hemi- 
spherical and fully emerged (Fig 1(b)). Between those 



regimes, the interface of the bath is slightly deformed by 
the bubble. The upper part of the bubble (above the 
meniscus) is approximated as a spherical cap This 
is merely due to a nearly perfect balance in hydrostatic 
pressure inside and outside the bubble (the weight of the 
film being negligible). The radius of this spherical cap 
is denoted R; it is the relevant parameter for the deter- 
mination of the lifetime. Since small bubbles are fully 
immersed, the upper part is not always visible (as in 
Fig 1(a) I . The radius is then estimated by the following 
pressure balance 



Aa/R = 2(7/ Rb + Pwgzs 



(1) 



where Rb is the radius of curvature at the bottom of 
the bubble, pw the water density, zb the depth of the 
bottom and a the interfacial tension. This balance is 
deduced from an integration of pressure variations along 
the dotted contour in FigO 




FIG. 3: PDF of the lifetime r for two given bubble volumes: 
(•) 6.07mm - (a) 25.1mm. Curves are guides for the eyes. 
The solid line is an extreme value distribution with a mean 
T,n = 57s and a standard deviation a = 12.6s; the dashed line 
is a WeibuU distribution with a scale parameter of 49.5s and 
an exponent of 1.73. 

The statistical distribution of the lifetimes r is shown 
for two given bubble radii: R ~ 6.07mm (0.5mL in 
volume) and R = 25.1mm (40mL in volume). A hun- 
dred bubbles have been studied for each size. In FigO 
the Probability Density Function (PDF) of the lifetime 
T has been plotted. For the 6.07mm bubble, the PDF 
is a right-shifted peak that is well-fitted by an extreme 
value distribution. Such a result has also been found in 
d, 3] . It suggests that the lifetime related to the peak 
is the determinist maximum lifetime Tm of the bubble. 
Lower lifetimes are due to random accidents. These acci- 
dents cannot occur after the determinist lifetime. For the 
25.1mm bubble, the distribution seems to behave like a 
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broad WeibuU distribution: no determinist trend can be 
emphasized. ExperimentaUy, the Hfetime of small bub- 
bles {R < 10mm) has been observed to roughly obey a 
right-shifted peak distribution, while the lifetime of large 
bubbles {R > 10mm) follows a broad and more or less 
uniform distribution. 




FIG. 4: Maximum lifetime Tm of the bubbles as a function of 
the radius R of the upper spherical cap. The line is related 
to the law described by Eq.Q. 



For each bubble size, about 10 measurements are made. 
When the lifetime r is predictable and obeys the first 
kind PDF, it makes sense to only consider the mean of 
the 5 greatest observed lifetimes in order to remove the 
premature accidental deaths. This averaged maximum 
lifetime r,„ is presented with respect to the spherical cap 
radius in Figd) For bubbles with a radius R < 10mm, 
T„i is roughly proportional to R. For larger bubbles, 
deviates from the linear law and the lifetime distribution 
becomes more spread. This lifetime variation with the 
bubble size is completely different from the one obtained 
with PDMS bubbles jl[ (the lifetime was inversely pro- 
portional to the radius). 

As mentioned in the introduction, the lifetime is re- 
lated to the thinning of the film located at the top of the 
bubble. The reasons for that decrease of the film thick- 
ness h are (i) the drainage by gravity and (ii) the evap- 
oration. Due to the high surfactant concentration, we 
can reasonably assume that the interfaces of the bubble 
are rigid ones. The situation is again different from the 
PDMS bubble [l| for which a plug flow is considered since 
the interface is not rigidified by surfactant molecules. At 
the top of the bubble, they are well approximated by a 
spherical cap of radius RzL h/2. 

The water film between both surfactant layers is thin 
enough to suppose that the flow obeys Poiseuille law. 
The tangential velocity U of this flow, averaged over the 
film thickness h, is given hy U = —{h^dep)/{l2fit^,R) 



where 9 is the azimuthal coordinate (starting from 
at the top), the dynamical viscosity of the water, 
and p the pressure inside the film. The pressure gradi- 
ent in the tangential direction is mainly due to gravity: 
dep — —pwgRsmO, where pw is the density of water. 

The drainage is particularly slow when the film is very 
thin. In these conditions, the evaporation process be- 
comes dominant compared to the drainage. Maxwell Q 
and Langmuir Q have independently proposed models 
of evaporation based on the diffusion of vapor molecules 
into the air. The radial mass fiux of vapor J is given 
by Pick's law as a function of the mass fraction of vapor 
rriv = Pv/ Pen where pv is the density of vapor and pa 
the density of air: J = —paDydniy/dr. In those models, 
the diffusion is supposed to be steady and the convection 
in the air is neglected. Practically, the first hypothe- 
sis is supported by the fact that the diffusion coefficient 
of water vapor in air is about Dy = 2.6 x lO^^m^/s. 
The relaxation time is about 4s for a centimetric bub- 
ble. The second hypothesis is roughly verified since 
in the experimental setup, the air velocity Va is very 
small near the bubble {Va ^ O.Olm/s). Ranz and 
Marshall [1^ have proposed a model to take convec- 
tion into account. According to them, the fiux is given 
by JR^{1 + QA2Re^/^Sc^/^)paDy[my{R) - m^i+oo)], 
where Re — VaR/va is the Reynolds number of the con- 
vection fiow, and Sc = Dyjva is the Schmidt number. 
Since at the interface, the vapor pressure p^ is equal 
to the saturation vapor pressure p^"', the difference in 
mass fraction can be expressed as a function of thermo- 
dynamical properties of vapor in air: TO„(i?)— mi,(-|-oo) — 
(1 - C){Mypl°-^)/{MaPa) where My and Ma are respec- 
tively the molar mass of vapor and air, pa is the par- 
tial pressure of the air and C is the relative humidity. 
From these considerations, the evaporation rate can be 
inferred: 



with 



pin Ma Pa ^ 



C)(l + 0.42i?ei/25ci/3 



(2) 



Since ~ O.QSlpa at 25C and ( = 0.54, the evapora- 
tion coefficient is about k ~ 2.8 x 10^^"m^/s. Due to con- 
vection, this coefficient is higher (about 0.5 x lO^^^m^/s) 
and it slightly depends on the Reynolds number (and thus 
on bubble radius R). 

The evolution of the film thickness is found by us- 
ing the continuity equation, and the following lubrication 
equation is obtained: 



dth^ 



smO \ R 



(3) 



As observed and explained in [ll|, the angular depen- 
dence is weak, in particular for 9 < tt/A: the film thick- 
ness is roughly constant in space over the upper part of 
the bubble. Moreover, the location of the minimum in 
thickness is obviously the top of the bubble. Therefore, 
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we can linearize Eq.(l3]) for positions near to = 0, lead- 
ing to: 

dh fc 

dt + 6^^R ^ R ^ ^ ' 

The maximum lifetime of a bubble is then given by 

R f'^° dh 

where /ip is the initial film thickness, the film thickness 
at rupture ^] and 0'^ = Q^^k/ Pwg- The [3 parameter is 
interpreted as the length scale at which evaporation and 
drainage are equally efficient. For h <^ /3, the film thins 
by evaporation much more than by drainage, while it is 
the opposite for 13. 

On a practical point of view, he ~ 10~''m while /3 « 6x 
10~^m. Therefore, /ic can be replaced by in Eq.®. In a 
similar way, when the initial thickness ho is much higher 
than /3, ho can be replaced by +00 without changing 
significantly the resulting lifetime r^. Therefore, one has 

~ ttR?! - (6) 

According to this expression, the lifetime r™ is roughly 
proportional to the radius of the spherical cap R (k 
slightly depends on R). This law is plotted in FigUl 
it is in good agreement with the experiments, at least 
for small bubble sizes. Note that no fitting parameter is 
needed. 

While ho ^ /3, the lifetime is predictable since it does 
not depend on initial and final thickness values. The de- 
viation to the theoretical law reflects the manner that we 
produce the bubble. To obtain a large bubble, some air 
is blown inside the liquid shell. The inflation thins the 
film; the initial film thickness /iq (after blowing) is thus 
lowered. When bubbles are large enough, ho < (] and 
it cannot be replaced by infinity in the lifetime compu- 
tation. The lifetime now depends on ho, and is smaller 
than its asymptotic value. Since ho is not controlled at 
all, T appears to be not determinist, as seen in the statis- 
tical analysis. The change in the PDF shape occurs for 
the same bubble size than the deviation from the linear 
scaling for t,„. 

In conclusion, an experimental work has been made in 
order to assess about the lifetime r of single soap bubbles 
created at the surface of a water bath. The maximum ob- 
served lifetime is shown to be roughly proportional to 
this radius for small bubbles. According to a statistical 
analysis, this maximum seems to be determinist and re- 
productible. Contrary to previous works, the evaporation 



seems to be the key physical process that determines the 
lifetime. The bubble is geometrically characterized by 
the radius R of the upper spherical cap that is formed 
when the bubble rises the surface of the bath. An analyt- 
ical model, based on drainage and evaporation of the wa- 
ter film under a small convection flow, has been proposed 
in order to explain this quasi-linear scaling. When the 
initial film thickness is large enough (typically for small 
bubbles), the lifetime does not depend anymore on ini- 
tial and final conditions: it is predictable. When bubbles 
are inflated too much, the initial film thickness becomes 
small enough to infer on the lifetime by reducing it. Since 
this thickness is not controlled, the lifetime looks unpre- 
dictable. When several bubbles are considered together 
(as in foams [l3|), Plateau borders probably play a sig- 
nificant role in the drainage process and the variation of 
the lifetime with the size is expected to be different. 
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